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The possibility of a node in the x dependence of the Sivers and Qiu-Sterman functions is discussed 
in light of its importance for the experimental check of the overall sign change of the Sivers effect 
between semi-inclusive DIS and the Drell-Yan process. An s-dependent version of the Ehrnsperger- 
Schafer-Greiner-Mankiewicz relation between the Qiu-Sterman function and a twist-3 part of g2 
is presented, which naturally suggests a node in the Qiu-Sterman function. This relation could 
be checked experimentally as well and could provide qualitative information on the gluonic field 
strength inside the proton. Satisfying the Burkardt sum rule by means of a node is briefly discussed 
and the importance of modelling the Sivers function including its full Wilson line is pointed out. 



PACS numbers: 13.88.+e,12.39.-x,13.85.Hd 

The Sivers [l| and Qiu-Sterman effects have been proposed as possible explanations of single transverse 
spin asymmetries Ajv observed in the process p -k X Q. In recent years it has become clear that these two 
effects are intimately related In this paper the x dependence of these effects is discussed, in particular the 

possibility of a node. 

The Sivers effect is described by a transverse momentum dependent parton distribution function (TMD) and 
generates azimuthal spin asymmetries for instance in the Drell-Yan (DY) process and in semi-inclusive DIS 
(SIDIS). The Sivers effect asymmetry in SIDIS has been clearly observed in the HERMES [| and COMPASS 
experiments. The Sivers function, here denoted by f^{x,k'^), describes the difference between the probability 
to find a quark with lightcone momentum fraction x and transverse momentum kx inside a hadron polarized 
transversely to its momentum direction and the one where the polarization points in the opposite direction. As 
the Sivers function describes a difference of probabilities it is not necessarily positive definite. In fact, the major 
interest in extracting the Sivers function from the DY process is that it is expected to have the opposite sign 
compared to the one extracted from SIDIS 3 ■ The gauge invariant definition of the Sivers function is in terms 
of a nonlocal operator involving a Wilson line: 

fl-^^\x,kl)efSTo.kT^ = yP.T. {P,ST\m^c{0,Ol^-m\P,ST}\^^^,, (1) 

where Cc denotes the Wilson line along contour C; St denotes the transverse spin vector; and F.T. denotes taking 
the Fourier transform, where ^~ and are the Fourier conjugate variables of xP~^ and /ct, respectively. The 
Sivers function is not uniquely defined, as it depends on the contour of the Wilson line, which in turn depends on 
the process considered. The Sivers function appearing in SIDIS contains a future pointing Wilson line, whereas 
in DY it is the same except past pointing, leading to the following overall sign relation [sl: 

/it H^T^rp) = —f-^r^ \x,krp). (2) 

This is a prediction of the TMD formalism that remains to be tested. A related sign test in W and Z production 
at RHIC has been put forward in Refs. [1, In more complicated processes that allow TMD factorization, 
Sivers functions with other Wilson lines can appear, which are not simply related by an overall sign to the Sivers 
function of SIDIS to which we will refer as "the" Sivers function from now on. 

It is important to emphasize that the above sign relation is about the overall sign and that the Sivers function 
itself need not be of fixed sign as a function of x. It can have one or more nodes in the x dep endence, even 



though its present extraction from SIDIS data in a restricted x range does not display a node Nevertheless, 
the possibility of a node should be kept in mind when comparing the extraction from SIDIS with the future one 
from DY. A node position is generally expected to be dependent, therefore, unless one compares the functions 
at the same x and values, the sign change test need not be conclusive. Moreover, such a node need not be 
at the same position for the different fiavors, possibly complicating the comparison further. Unless one makes 
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sure that nodes do not play a role in the comparison, one may wrongly jump to the conclusion that the TMD 
formalism is flawed in some way if the overall sign change between SIDIS and DY is not confirmed in experiment. 
Motivated by the importance to know whether the Sivers function has a node, we investigate if there are any 
other indications in favor or against such a node. 

Since we are interested in the x dependence here, we will not address the transverse momentum dependence 
(for a discussion of possible nodes in the kx dependence cf. and restrict to the first transverse moment of 
the Sivers function, i.e. the Sivers function weighted with the transverse momentum squared: 

This quantity is of interest because of its direct relation to the twist-3 Qiu-Sterman (QS) function T{x, St) 0^ 

T{x,St)^i^ I '^e^^^P,S\mV^j dr^F+'^irjn.) ^(An_)|P,5), (4) 

where Fq, = eTpaS^fi- / {2iM P^) (like in [2] we take 5^ = 1) and e^" = e"^'"^n_|_ctn_^. The relation between the 
first moment of the Sivers function (of SIDIS) and the QS function is a direct proportionality [4]: 

ft'^'\x)=-^T{x,ST). (5) 

The a:-dependence of the two functions is therefore the same, apart from an overall proportionality constant. 
From now on we will absorb the coupling constant g into the definition of T, but it is displayed here explicitly 
since it determines the relative sign in front of it (for a discussion of this issue cf. J^] ) . 

The above definition of the QS function is given in the ~ lightcone gauge for simplicity and contains 
the operator J drj F~^°'(rin-), which renders it intrinsically nonlocal along the lightcone, even in the = 
gauge. Below we are going to discuss two different assumptions about this lightcone integral over the gluonic field 
strength, one by Qiu and Sterman and one by Ehrnsperger, Schiifer, Greiner and Mankiewicz (ESGM) [l6| . 

Qiu and Sterman considered the following parametrization [2'|: 

T«(x,5T) = «,A/f(x), (6) 

where fi{x) is the ordinary unpolarized parton distribution function and q denotes the quark flavor. This 
parametrization follows when the QS matrix element is viewed as yielding the average value of J drj F^°'{'i]n^) 
inside the unpolarized proton and hence is expected to be simply a number times the unpolarized distribution 
function: 

/i(^) - ^/ ^e''^(P|?^(0)^_^(An_)|P). (7) 

The function /i(x) is a probability distribution and hence of definite sign. Using this parametrization, T{x, St) 
can have different signs for different flavors, but cannot exhibit any node. In order to roughly describe the single 
transverse spin asymmetries (SSA) experimentally measured in the process p ^ tt X a.t y/s ~ 20 GeV [3| , 
the following choices were made for the parameters: Ku = +1 = —Kd, Ks = 0, yielding A ~ 100 MeV ^. This 
parametrization was subsequently used to predict SSA in tt production at y/s = 200 GeV [13] and in Drell-Yan 
[ij, [13 • No conclusive evidence in favor of the above parameterization has been obtained yet t houg h. 

Another view on the lightcone integral over the gluonic field strength is taken by ESGM. In Ref. [1^ it is assumed 
that the field strength is only significantly contributing inside the proton and that the following approximation 
holds: 

dr; F+"(77n_) « F+"(0) X 2cMi?o, (8) 

where Rq is the proton "lightcone" radius in the rest frame of the proton (taking into account that it is being 
probed by a highly relativistic probe, which sees t he p roton as Lorentz contracted) and V = J dj] = 2cMRo, where 
c is argued to be a constant between 1/3 and 1 jl6l. fT7l|. ESGM considered the field strength at 77 = 0, because 
that corresponds to the position of the quark fields upon integration of T(x, St) over x, such that a local operator 
is obtained. A priori it is not known whether this approximation is legitimate, except that it may simply be true 
numerically for some configurations, but of course there is no way to select such configurations. Nevertheless, 



3 



it leads to an interesting result, which is the ESGM relation between the lowest (zeroth) Mellin moment of the 
QS function and the second moment of the twist-3 part of the distribution function 52, which according to 1^ is: 



1 pI 

2d / ^2 



T{x,St) dx ^ -l2cM Ro X ff2(a;)ltwist-3 ^2;. (9) 

Jo 

It should be emphasized that this is not an exact relation. But if ESGM's approximation is fine, then one 
can relate the magnitude of the above mentioned SSA to g2] this would be very useful, since it is known that 

/ x'^ 52(2;) I twist-3 is very small. 

The structure function 172, which in the parton model is directly related to the distribution function g| via 
g2{x) = ^ q ^q92' with 6^ the quark charge squared in units of the electron charge, has been measured by the 
E155 experiment at SLAG. Also the twist-3 part of g2 was extracted, yielding a value for its second moment d2 
which is defined as 

da = 3 / x^ 52 (a;) I twist-3 dx (10) 

The E155 experiment obtained d2 — 0.0032 ± 0.0017 for the proton [l^, when taking into account all available 
SLAG data. Assuming there are no large cancellations among the quark fiavors, one concludes that also c?2 and 
(^2 are both very small. This conclusion is supported by lattice QGD evaluations [T9l |. 

Gombining the ESGM relation and the above parametrization ([6]) of the QS function would suggest a very small 
SSA in p TV X contrary to observations. Given the very small size of c?2, one would conclude that either A is 
much smaller than expected from the SSA data or Rq has to be unnaturally large 1 fm). Failure of the ESGM 
relation is one possibility and may indicate some qualitative features of the lightcone integral of the gluonic field 
strength inside the proton, such as that it is not slowly varying inside the proton or is significantly contributing 
outside the proton too. But one could also question the validity of parametrization because one way to allow 
for large SSA and simultaneously accomodate small ^2 through the ESGM relation is to consider the option that 
T{x, St) changes sign as a function of x, having large absolute value in certain x regions, but having a small 
integral. This is the view we will explore here, leaving aside the idea that T{x, St) is proportional to fi{x). 

The approximation ([5]) can be extended in various ways to obtain an x-dependent version of the ESGM relation. 
For instance, assuming that the gluonic field strength is a slowly varying function (or constant even), within the 
proton, leads to the approximation 

drj F+'^{r]n^) « F+"{r]an^) x 2cMRo, (11) 

for some (or any) rjo within the proton lightcone radius Rq. To reduce to the ESGM relation upon integration 
over X, one has to consider rjo at the position of either one of the quark fields. This is automatically satisfied if 
the field strength is taken to be constant within the proton. It should be mentioned that if one views the above 
approximation as giving the average field strength times the integration region, that in that case the estimate of 
c of Refs. [l^, [l3| rnay be far off if the gluonic field is in fact heavily fiuctuating. This can be experimentally 
investigated. 

With the approximation (|lip for rjo = A, the following unintegrated relation can be obtained^: 

T'^ix, St) = -2cM'^Ro x^ g^{x). (12) 

This is an unintegrated version of the ESGM relation, which holds for each flavor separately. The twist-3 distri- 
bution function gT{x) is the quark-gluon-quark correlation part of gT{x) = gi{x) + g2{x) split off by means of 

the equations of motion (for m = 0): gT{x) — gT^x) — g[}^{x)/x (cf. e.g. Refs. (2l| - [23| V It can for instance be 
measured in the Drell-Yan process in double spin asymmetries Alt of longitudinally polarized hadrons colliding 
with transversely polarized hadrons [l^, [2l| or using SIDIS [1^ . 
Upon taking the lowest Mellin moment of Eq. ([T^. one arrives at: 

J T'^ix, St) dx = -cM^Ro [d^ + 4] , (13) 



To obtain this result the lightcone integral is implemented as centered around zero, such that in the notation of _2(J the correlator 
#^(a;,3/) is symmetric under the interchange of x and y and ^'^{x,y) antisymmetric and therefore GAi^^y) = 0. 
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where 

dl = 3 x"^ gl{x)\^^i^^_^dx = 2 x'^g^x)dx. (14) 
Jo Jo 

Summing over quark flavors, one obtains a relation in terms of the twist-3 part of the structure function g2'. 

Ve, / T'^{x,St) dx = -6cM^Ro (a;) | twist-3 (15) 

g J-i Jo 

which apart from the sum over flavors and the quark charge squared factor is a factor of 2 different from Eq. ^ . 
Given the uncertainty in the proportionality constant c this factor is not of importance, but we do note that the 
flavor dependence of the ESGM relation was not addressed properly in Ref. pjj . 

The question here is whether grix) has a sign change as function of x, since the unintegrated "ESGM" relation 
implies similar behavior for T{x,St)- The bag model suggests that 52(2;)! twist-3 ^ ^^S^ changing function 
of X. However, ^2 (a; )| twist-3 ^^"^ drix) correspond to different operator matrix elements, even though the second 
moments are directly related through Eq. (HH). Instead, one can look at the first moment of grix), which in the 
= gauge is given by: 

'^xgUx)dx = Re (p, 51^^(0)7^755^^(0) V(0)|P, 5)^^^^ = 0, (16) 

where the vanishing of this particular matrix element is shown in Ref. [25l | on the basis of Lorentz invariance, 
analogous to the derivation of the Burkhardt-Cottingham sum rule <?2(^) ~ 0- Since the first moment of 
grix) deals with a local operator that appears in the Operator Product Expansion (OPE) for charged currents, one 
can more specifically conclude that also Jp xg!^{x) dx corresponds to a local operator matrix element (obtained 

through Taylor expansioii) that vanishes due to Lorentz invariance, in contrast to gKx) da; which does not 
appear in the local OPE 0, [i?!- In other words, the vanishing integral in Eq. (fTO]) is not due to a cancellation 
among the quark {x > 0) and antiquark {x < 0) contributions. 

The vanishing of x g^{x) dx implies that g^{x) has a node and through the unintegrated "ESGM" relation 
also the QS function would have a node. In that case, J T{x, St) dx can be much smaller than the maximum 
value of T{x, St), such that a sign change of T{x, St) can accomodate both small d2 and large SSA in a limited x 
range in a natural way. The unintegrated relation is useful in this respect, since asymmetries can be measured as 
function of a;, therefore, the relation can be checked in experiment. The option that the QS function has a node 
should then be kept in mind and also that the node can change position as a function of Q^, which is relevant 
when comparing different experiments. 

What speaks against a node is that most model calculations of the Sivers function (28l - [33 | do not show a node 
(ignoring some small bag model artifacts), except for those of Refs. [H, [11] which obtain down quark Sivers 
functions with a node, albeit very different ones, and for Ref. [3^ which obtains an up quark Sivers function with 
a node. These model calculations of the Sivers functions all consider the gauge link to lowest nontrivial order in 
the coupling constant, in other words, the first order expansion of the Wilson line. It is unclear what the size of 
the higher order corrections is and whether these could change the sign in a particular x region. 

A common feature of the model results is that the up and down quark Sivers functions are opposite in sign. 
This is also expected from the large Nc limit, which leads to f^{x, k"^) = —f^{x, k^) + 0{N~^) [s^, [s^, and 
on the basis of the signs of the quark anomalous magnetic moments through an integral relation of the Sivers 
function and GPDs 40, 41]. This suggests that if a node occurs, it is present in both up and down quark Sivers 
functions. 

SIDIS experiments off proton and deuteron targets in the valence region seem to indicate that up and down 
Sivers functions are indeed opposite in sign and moreover similar in magnitude. In some models the magnitude 
of the up quark Sivers function was found to be much larger than that of the down quark (sol . |35| , but in more 
recent model calculations they are found to be of comparable magnitude [Sll [ssl Is^ . [361 [37| . If indeed of opposite 
sign, but similar in magnitude, the so-called Burkardt sum rule |42l | 



a=q,g 



ft^^'^^x) dx = 0, (17) 



5 



can be satisfied by a cancellation among contributions of the valence quarks, not requiring large contributions 
(and accompanying large cancellations) from non-valence contributions. This would be one possibility. See Ref. 
[43.] for a check of the Burkardt sum rule in a model calculation that includes the gluon Sivers function. A second 
possibility to satisfy the Burkardt sum rule is that the x integral of the Sivers moment is zero for each flavor 
separately. The latter is not in contradiction with the measured Sivers asymmetries in SIDIS, since those do not 
provide full information on the integral over x of the Sivers moment. A future observation of a node in the Sivers 
function could point to this second possibility, as does a small gluon Sivers function. The latter can be seen by 
rewriting the Burkardt sum rule in the form: 

V/ T''{x,St) dx = - f T^ix,ST)dx, (18) 
g J-i J-i 

and by comparing it to the integrated ESGM relation of Eq. ((T5|) . where the l.h.s. of the latter includes an 
additional quark charge squared factor. If the gluon QS (or Sivers) function turns out to be small, then combined 

with the fact that ^2 is small, it implies that both J2q I^i St) dx and Yl,q lli St) dx are small, 

disfavoring large cancellations among the flavors, but rather suggesting the x integral of the QS function to be 
small for each flavor separately. In this way a small gluon QS/Sivers function and a node in the quark Sivers 
functions would go hand in hand. 

In conclusion, the x-dependent version of the ESGM relation presented here implies that the QS function changes 
sign as a function of x, but it is based on some assumptions that may not hold. It would be very interesting 
to test this relation in experiment, even outside the region of a possible node, since the proportionality relation 
should hold for all x values. Confirmation of or deviations from the relation could teach us about the qualitative 
features of the giuonic field strength inside the proton. Despite this uncertainty a node should be considered as 
a serious possibility, since twist-3 functions are not probability densities, therefore need not be positive definite. 
It would be just as natural for the QS function to have a node as it is for the twist-3 part oi g2- A sign change 
as function of x also offers another way to satisfy the Burkardt sum rule, without requiring large cancellations 
among the different flavors. This sum rule together with the integrated ESGM relation obtained here suggest that 
a small gluon QS/Sivers function disfavors such large cancellations. A node is found in a few model calculations 
for the Sivers function, but not in most, and not for up and down quarks simultaneously, contrary to expectation 
from the large Nc limit. Extending these model calculations of the Sivers function to include its full Wilson line, 
or equivalently performing model calculations of the Qiu-Sterman function, would therefore be very interesting. 
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